Giant Magnetoresistance in Multilayers with Noncollinear Magnetizations 
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We study the dependence of perpendicular-current magnetoresistance in magnetic multilayers on 
the angle between the magnetizations of the layers. This dependence varies with the thickness of one 
of the layers, and is different for multilayers with two and three magnetic layers. We derive a system 
of equations representing an extension of the two-current series resistor model, and show that the 
angular dependence of magnetoresistance gives information about the noncollinear spin-transport 
in ferromagnets. 
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The discoveries of giant magnetoresistance (GMR) Q 
and spin-transfer Q in ferromagnetic metallic multilayers 
greatly contributed to our understanding of the relation 
between magnetism, charge, and spin transport, and lead 
to important applications in memory devices and sen- 
sors. The spin-torque theory of spin-transfer relies on 
the absorption of the transverse spin-current at the mag- 
netic interfaces in multilayers with noncollinear magneti- 
zations, due to the averaging of spin-dependent electron 
reflection at the interfaces, and spin- pr ecession inside the 
ferromagnets I1I1I1I1I1I1I1000 Similarly the 
disappearance of the transverse spin-current inside ferro- 
magnets is predicted to affect the angular dependence of 
perpendicular-current (CPP) GMR in multilayers with 
noncollinear magnetizations (AGMR). Thus, AGMR is 
an important effect complimentary to the spin-torque. 

Theories of AGM R q ualitatively agree with the avail- 
able data 0,0,0,0], but quantitative agreement has 
not yet been achieved. Due to the lack of systematic stud- 
ies, it is also impossible to verify the predicted tendencies 
for the variation of AGMR with the multilayer parame- 
ters. Here, we present a systematic study of AGMR in 
multilayers with two and three magnetic layers, in which 
we varied the thickness of one of the layers. The depen- 
dence of AGMR on this thickness is different in the two 
studied structures. Our analysis shows that magnetic in- 
terfaces, transverse to magnetizations spin-currentsw in 
ferromagnets, and in some cases sample leads give con- 
tributions to AGMR. These findings are important for 
the understanding of spin-transport in ferromagnets and 
theories of spin-torque. 

Our sample fabrication and measurement techniques 
were described elsewhere (ljj . The structure of sample 
type A was Nb(150)Cu(20)FeMn(8)Py(6)Cu(10)Py(1.5- 
12)Cu(20)Nb(150), Py=Permalloy=Ni 84 Fei 6 (Fig. [TJ. 
We label these samples A1.5, A3, A6, A12 by the 
thickness tp y of the top Py layer. All thicknesses 
are in nanometers. Samples B 1.5- 12 had struc- 
ture Nb(150)Cu(20)FeMn(5-12)Py(6)Cu(10)Py(1.5- 
12)Cu(10)Py(6)FeMn(5-12)Cu(20)Nb(150), and were 
labeled by the thickness of the middle Py layer. The 
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FIG. 1: (a,b) Schematics of samples A (a) and B (b), with 
layer compositions as labeled (thicknesses not to scale). Dif- 
ferent sections of samples are labeled 1-3, as described in the 
text, (c) Schematic of the model incorporating finite trans- 
verse spin-current decay length. 



bottom Py(6) in samples A, and the outer Py(6) layers 
in samples B were exchange-biased by annealing in 
magnetic field H « 30 Oe at 170°C. At least 2 samples 
of each type were measured with similar results. 

Sample resistances R were measured with a SQUID 
voltmeter at 4.2 K in a CPP geometry. Fig. Ota) shows 
an example of R vs. H curve for sample A12, measured 
with H > along the pinning direction of the bottom 
Py layer. At large H > 0, magnetizations of both Py 
layers are parallel (P) to each other and H, giving low 
resistance Rp. As H decreased to small H < 0, the free 
layer switched to give antiparallel (AP) state with high 
resistance Rap- The free layer coercive field H c (half the 
width of its hysteresis) decreased from sw 30 Oe in sam- 
ples A1.5 to 10 Oe in A12. The decrease of resistance to 
Rp at significantly larger H < is due to switching of 
the pinned Py layer. We define the bias field Hb as H 
necessary to achieve R = (Rp + Rap)/ '2. In most sam- 
ples, \Hb\ ~ 800 Oe. Because H c <c Hb, we were able to 
rotate the free Py layer's magnetization, not significantly 
affecting the pinned one. 

Fig.^b) shows examples of the AGMR measurements 
for samples A1.5 and A12, performed by rotating a fixed 
H = 25 — 100 Oe in the plane of the films. At angle 
8 = 0, H is in the pinning direction of the bottom Py. 
There was no significant dependence of data on 25 < H < 
100 Oe in samples A3-A12. In samples Al.5, H — 25 Oe 



2 



14 



a 
en 



13 



(a) 

7 / 




(b) rfR 

r V 




















-1.0 -0.5 0.0 0.5 


100 200 


H (kOe) 


8 (deg.) 



FIG. 2: (a) i? vs. # in sample A12. (b) AGMR at H = 25 Oe 
(dots) and H = 100 Oe (squares) in samples A1.5 (bottom) 
and A12 (top). Solid lines are fits of the 100 Oe data with 
X = 7.7±0.6 and % = 1.96±0.05, respectively, as described in 
the text. For clarity, not all the measured points are shown, 
and the A1.5 data are offset by — 0.2nf2. 



was insufficient to completely reorient the Py(1.5) layer 
due to its higher coercivity, but at 50 < H < 100 Oe 
data were independent of H. Thus, we conclude that: a) 
our H < 100 Oe <C Hb does not significantly affect the 
magnetization of the pinned layer, b) Except for samples 
A1.5, H — 25 Oe was sufficient for monodomain rotation 
of the free layer. 

Solid lines in Fig.^b) are fits of the 100 Oe data with 



R(9) = Rp + AR 



1 -cos 2 (6»/2) 
l + XCOs 2 (0/2) : 



(1) 



proposed by Giacomoni et al. 15 L and later derived for 

symmetric spin- valves II II II Oil II H El in E q □ 

AR = Rap — Rp 7 and x is a fitting parameter. Eq. 
gives good fit for all samples except for A1.5. In all three 
samples A1.5, the best fit did not completely reproduce 
the nearly constant data for —90° < 9 < 90°. Because of 
the finite 9 = curvature of the fit, it was below the data 
at 9 « 0, and slightly above it around 9 = 100°. Data for 
Al.5 in Fig. I^b) show a weak rise around 9 = 0, similar 
to the uncertainty of the measurements. 

Both R(H) and AGMR data for samples B were quali- 
tatively similar to those in Fig. [5] Figs.|Hta,b) summarize 
the AAR measurements (where A is the sample area) for 
samples A and B, and Figs. EUc,d) show the x values, 
extracted by fitting the AGMR data with Eq. |T| To re- 
duce the uncertainty, x obtained from 100 Oe and 50 Oe 
data were averaged for each point. At large tp y , sample 
B is expected to behave similarly to two samples A con- 
nected in-series. Consistently, Fig. shows that AAR 
for samples B12 is double that for samples A12, while x s 
are close. At smaller tp y , AAR in samples A decreases 
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FIG. 3: (a) Measured (symbols) and calculated (line) AAR 
vs. tp y for samples A. (b) Same as (a), for samples B. (c) 
X from the fits of measured (symbols) and calculated (lines) 
AGMR with Eq.^for samples A. Solid line is calculated with 
correction for a finite l t = 0.8 nm, dashed line - with l t = 0. 
Inset: calculated R(9) for Al.5, as described in the text, (d) 
Same as (c), for samples B. All the lines are B-spline fits of 
the calculated points for tp y = 1.5, 3, 6, 12 nm. 



faster than half of that in samples B. In samples A, x 
increases, while in samples B it decreases at smaller tp v . 

Several models have treated AGMR 0, H S 0, S H 
EIHII3- The y predict the form Eq.|T]for symmetric 
spin-valves, but overestimate x typi ca Uy by a factor of 
two. Some theories predict a negligible length scale It 
for the decay of the transverse spin-current in Py 0, . 
Others correlate It with the magnetic length « 4 nm in 

Py 13 

We develop an extension of the two current series resis- 
tor model [13 (2CSRM) to AGMR, which contains much 
of the same physics as these other models, but allows us 
to gain qualitative insight into our results, analyze the ef- 
fects of finite It, and include the effect of the Nb sample 
leads, which we shall see give an important contribution 
to the measured dependencies. 

In 2CSRM, one separately considers two spin-current 
channels, same across the whole sample. For noncollinear 
magnetizations, two common spin-channels across the 
multilayer generally do not exist, so one needs to con- 
sider how the spin-channels are transformed across the 
multilayer. In samples A, we consider three parts of the 
multilayer (labeled 1-3 in Fig. ^a)), separated by the 
Py/Cu interfaces 12 and 23. The outer limits of regions 
1 and 3 are determined by the spin-diffusion, i.e. the 
GMR-active part of the multilayer. This additional con- 
straint augments 2CSRM and is essential for the follow- 
ing analysis. At each point of the multilayer, we define 



3 



a matrix current / = 



which gives charge 



pi pa 

j21 j22 

current / c = TrI and projection of spin current on an 
arbitrary axis m, I s (m) = TiTVQ^j rnidil) / (2e) . Matrix 
product of / and the Pauli matrices er^ is implied here. 
We assume that, due to averaging of spin-precession in 
the ferromagnets, the component of spin-current trans- 
verse to the Py magnetizations vanishes in regions 1 and 
3 Later, we will discuss and modify this assump- 

tion. In region 1, the matrix current I\ is then diagonal 
in the frame set by the magnetization M\ of the bottom 
Py layer, I\ = diagl\. Similarly, in region 3, I' 3 = diagl^. 
We use primed symbols for the frame set by the orienta- 
tion of the top Py magnetization M 2 . 

If spin-flip scattering within regions 1-3 is neglected, / 
is conserved separately across each region. However, only 
I c and the spin-current projections on the corresponding 
Py magnetizations are conserved at the interfaces 12 and 
23, therefore I\ = diagl 2 , I' 3 — diagl'%. Here V = 



U+IU, U 



is the spin rotation 



cos(0/2) -sin(0/2) 
sin(0/2) cos(0/2) 
matrix by angle between Mi and M2. 

We describe the local electron distributions by 2 x 2 
spinor distribution matrices 0, Ej. Their diagonal 
elements in an arbitrary reference frame are given by 
the spin-up and spin-down electron densities, related 
to this frame. We neglect scattering in the Cu spacer 
between the Py layers, i.e. assume a position indepen- 

Finally, we introduce 
1-/3! 



dent distribution in 
matrix resistances 



R 



2R* 12 



2R^U 



region 2. 

Ri = 
1 



2R{ 
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and R23 = URi2U + (for identical 
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interfaces 12 and 23), which connect martrix currents 
across regions 1 and 3, and interfaces 12 and 23, to the 
corresponding variations of the electron distributions. 
R* , 7, and (3 are standard GMR notations |l7|. 



The total voltage across the multilayer V — 



is 



the sum of electron distribution variations across regions 
1 and 3, and interfaces 12 and 23, 

V = Rih + (-R12 + #23)^2 + Rsh = Ridiag{I 2 ) + 

(R12 + R.23% + R 3 Udiag{U + i 2 U)U + . (2) 

Eq. [21 connects four unknown components of I 2 to the 
voltage V across the mutlilayer. The diagonal compo- 
nents of Eq. reduce to 2CSRM in the collinear limit. 
Once Eq. is solved for I 2 , the resistance of the mul- 
tilayer is given by R(6) = V/Trlz. We note that while 
Eq.[5]is written in the frame of M\ , it can be transformed 
to an arbitrary reference frame. 

For a symmetric bilayer, R\ = R£, 0i = 0a, Eq. [21 is 
diagonal in the frame rotated halfway between M\ and 



M.2- In this special case 
R(0) = 2R{ + 2R* 12 



2cos 2 (0/2)(R* 1 p 1 +R* 12 j) 2 



R 



(3) 



12 



Rl cos 2 (0/2) 

This expression has the same form as Eq. ^ with x = 
R*/R* 2 . As we noted above, Rl and 0\ are determined 
by the condition that only the GMR-active part of the 
multilayer is included. In particular, in samples A12 
tp y > l a f ss 5.5 nm, the spin-diffusion length in Py |l8j| . 
so only the Z s /-thick parts of the Py layers should be in- 
cluded in regions 1 and 3. Thus, samples A12 are approx- 
imately symmetric, and can be described by Eq. [31 Sim- 
ilar arguments hold for samples B12, viewed as two sam- 
ples A connected in series. We use the parameters estab- 
lished in CPP-GMR measurements [II], AR\ = l s fP* Py = 
1.4 mm 2 , Pi = l3 Py = 0.7, AR* 12 = 0.5 ffhn 2 , 7 = 0.7. 
We obtain x — R*/R*2 = 2-8, larger than the measured 
X = 2.0 for A12, and x = 1-6 for B12. 

Overestimation of x is a general tendency of the 
AGMR theories. The physical transparency of Eq. [3| al- 
lows us to identify the possible origins of this discrep- 
ancy, and appropriately correct our analysis. We as- 
sumed above that transverse spin-current and the nondi- 
agonal electron distribution components vanish in Py ar- 
bitrarily close to the Py/Cu interface. However, a sharp 
Py/Cu interface may be an unjustified idealization. If 
the onset of the bulk ferromagnetic Py properties occurs 
over a finite thickness, where Py and Cu are alloyed, 
it is also reasonable to expect that the transverse spin- 
current decays over a finite length Z t , nominally inside Py. 
Moreover, some theories predict a finite l t even when an 
ideally sharp interface is assumed ■ Regardless of the 
physical origin, we can phenomenologically include a fi- 
nite It into our model by expanding the interfaces 12 and 
23 into finite /j-thick regions of Py (Fig. [He)), where 
spin-current is noncollinear to the magnetization. This 
correction decreases Rl, and increases R\ 2 , thus decreas- 
ing x- U ~ 0.8 nm gives a good agreement of calculated 
X with the data for samples A12 and B12. 

AAR is not affected by the finite It, irrelevant for the 
collinear transport. For tp y comparable to It, the spin- 
torque should decrease due to incomplete transverse spin- 
transfer between electrons and magnetization. We note 
that in the published studies of spin-transfer (mostly 
with Co) the ferromagnet thicknesses were larger than 
h UM- It in Co is likely even smaller than in Py due 
to its larger exchange splitting 7]. The circuit the ory 
of spin-polarized transport and spin-transfer 

El HUBS 

uses a mixing conductance parameter gn, characterizing 
spin-dependent scattering at the interfaces. Related to 
our model, 2g^ = (R! )~ 1 , thus it depends on l t . 

Giacomoni et al. 15] obtained x ~ 1-2 for thick 
Py/Cu/Py spin- valves similar to our samples A12 or B12. 
Their smaller x ma y be due to stronger alloying of the 
Py/Cu interfaces (giving larger l t ), caused by higher de- 
position rates and longer annealing time during pinning. 
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We modeled the \ dependence on tp y by solving the 
general form of Eq.[5] The solution for asymmetric multi- 
layers is similar to Eq.^ but more complicated, and the 
denominator contains additional cos 4 (0/2) terms. We 
therefore give only numerical results for specific cases. 
One has tp y < l s f in samples A1.5 and A3, so the GMR- 
active top part of the multilayer (region 3) must now in- 
clude the entire top Py layer and the Py/Cu and Cu/Nb 
interfaces. Studies of Nb/Py and Nb/Cu/Py multilayers 
yield a large interface resistance ARp y /Nb ~ 3 film 2 [18| . 
both with and without a Cu spacer between Py and Nb. 

For samples A1.5, we added the full value of Rp y /Nb to 
i?3, neglecting electron spin-flipping in Py(1.5). For A3, 
we added 0.5i?p a /jvb to approximately account for spin- 
flipping in Py(3), reducing the contribution of Py/Nb 
interface to GMR. Our model gives the same results for 
samples A12 and A6, but the difference in data is also 
small. The calculated AAR (Fig. |3a), solid line) agrees 
well with the data. For samples A1.5, the calculated 
R{9) significantly deviates from the form Eq. Inset in 
Fig. |3{c) shows that it has maxima both at 6 = and 
9 = 180°. A similar behavior is predicted for asymmet- 
ric spin- valves by the circuit theory . Our calculation 
exaggerates a weak rise of data at 9 ss (Fig. Gib)), but 
captures the overall experimental behavior. The quanti- 
tative discrepancy may be due to the neglected electron 
spin- flipping at the interfaces. In Fig. E^c), solid line 
shows the calculated \i defined as the best fit with Eq.Q] 
to the calculated R(9), using l t = 0.8 nm. The calcula- 
tion using It — (dashed line) gives a significantly worse 
agreement with data. 

Samples B are symmetric with respect to the center 
of the middle Py layer. Therefore, current reversal does 
not change the electron distribution at that point for any 
magnetic orientation of the middle Py layer. Since the 
properly offset electron distribution is proportional to the 
current, we conclude that in the center of the middle Py 
layer the electron distribution is spin- independent. The 
model developed above for two magnetic layers can now 
be adopted to samples B with tp y < 2l s f, if we take half 
of the middle Py layer as region 3, as shown in Fig.^b). 
The top half of the sample simply doubles the resistance 
obtained from Eq. not affecting \- The results for 
AAR(tp y ) and x(i_pj,), with l t = 0.8 nm, are shown 
with solid lines in Figs. EHb,d). The deviations from the 
form Eq.^were negligible for all samples B1.5-B12. Our 
model overestimates AAR, but gives reasonable results 
for x(tpy)- Calculation assuming l t = (dashed line in 
Fig. Etc)) gives a worse agreement with data. 

Qualitatively, our results for x(tpy) m both samples A 
and B can be understood with Eq. derived for sym- 
metric multilayers. In samples A, the activation of the 
highly resistive Py/Nb interface at smaller tp y is equiv- 
alent to an increase of i?* in Eq. |3 giving larger \- I R 
samples B, smaller tp y is equivalent to reduced R\ in the 



symmetric case, and thus smaller x- 

In summary, we showed that the variation of GMR 
with angle between the magnetic layers (AGMR) depends 
on the thickness of one of the magnetic layers. The de- 
pendence is different in samples with two and three mag- 
netic layers. To analyze the data, we developed an ex- 
tension of the two current resistor model to multilayers 
with noncollinear magnetizations. Our analysis leads to 
the following conclusions: i) The deviation of AGMR 
from sinusoidal behavior is approximately given by the 
ratio of two quantities: a) the resistance of the GMR- 
active part of the multilayer excluding the noncollinear 
ferromagnetic interfaces, b) the resistance of these in- 
terfaces. The magnitude of this effect is not directly 
related to magnetic anisotropics and AAR; ii) AGMR 
can be nonmonotonic in asymmetric spin-valves; iii) the 
transverse spin current penetration length ti into ferro- 
magnet can be extracted from AGMR. From our model, 
It ~ 0.8 nm for Py. l t is an important parameter for 
the models of noncollinear spin transport in ferromag- 
nets and spin-torque. The spin-torque should be reduced 
if the ferromagnet thickness is close to It- 
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